Effect of zero point phase fluctuations on Josephson tunneling 
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In the presence of phase fluctuations the dc Josephson effect is modified and the supercurrent at 
zero voltage is replaced by a peak at small but finite voltages. It is shown that at zero temperature 
this peak is determined by two complementary expansions of finite radius of convergence. The 
leading order expressions are related to results known from the regimes of Coulomb blockade and 
of macroscopic quantum tunneling. The peak positions and the suppression of the critical current 
by quantum fluctuations are discussed. 
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The dc Josephson effect allows a Cooper pair current 
to flow through a superconducting tunnel junction in the 
absence of an external voltage. The current is determined 
by the difference ip of the condensate phases on the two 
sides of the junction through I = Ic sin((/5) and is lim- 
ited by the critical current Ic- While this feature in the 
current- voltage characteristic has zero weight, it acquires 
a finite width due to either thermal or quantum fluctua- 
tions of the phase difference. Here, we analyze the role of 
quantum fluctuations and concentrate on the Josephson 
peak at zero temperature. 

The finite capacitance C of a Josephson junction is a 
source of fluctuations of the phase difference (p and thus 
of the broadening of the Josephson peak since the charge 
Q on the capacitance is the conjugate variable to (p [yj. 
At finite voltage V a current can only flow if the tun- 
neling Cooper pairs can loose their excess energy 2eV. 
This energy can be transferred to the degrees of freedom 
present in the electromagnetic environment of the junc- 
tion described by the impedance of the circuit. In order 
to observe a peak in the current as a function of the volt- 
age, the Josephson junction needs to be voltage-biased. 
This has become possible experimentally only very re- 
cently Since we are interested in the behavior close 
to the ideal Josephson peak, typical voltages V are much 
smaller than the gap voltage A/2e and quasiparticle ex- 
citations can be neglected. 

A minimal circuit displaying quantum fluctuations of 
the phase is shown in fig. |l| and may be described in terms 
of the Hamiltonian [|j 
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The first term corresponds to the charging energy and 
introduces an energy scale Ec = 2e^/C. The second 
term describes the tunneling of Cooper pairs through 
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FIG. 1. Josephson junction characterized by the Joseph- 
son energy Ej and capacitance C coupled to an ideal voltage 
source V via an external impedance. 

the junction. The Josephson coupling energy Ej is re- 
lated to the critical current by Ic — 2eEj/h. The third 
term describes the coupling of the junction to an exter- 
nal impedance modeled by a set of LC-circuits and also 
takes into account an applied voltage V. 

Summing the perturbative expansion in the Joseph- 
son coupling to all orders, one finds for the equilibrium 
Cooper pair current fl 



2n 



{C:»d 



k=l 



X 



dti 



dt2n-l GXp(r) 



(2) 



/O ^0 

where the exponent in the integrand is given by 
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The sums over Cfc and rjk run over the values ±1 with the 
constraint X]fc"o^ ^fc ~ 

This result depends on the phase autocorrelation func- 
tion in the absence of tunneling |^ which at zero tem- 
perature is given by 
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It is entirely determined by the total impedance ZtiLo) = 
[iu)C + 1/Z{iu)]~^ seen by the junction. In the following, 
we assume a purely ohmic external resistance R thereby 
neglecting features in the impedance, like resonances in a 
transmission line, which might lead to additional struc- 
ture in the J-F-curvc ||J^. One then has 
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where p — R/Rq with the resistance quantum Rq — 
h/Ae^. The total impedance is cut off at a frequency 
ijjji = 1/RC due to the junction capacitance. For suffi- 
ciently long times the correlation function is given by 



J(t) = -2p ln(wfl|t|) +7 + i-sign(t) 



(6) 



where 7 = 0.5772 ... is the Euler constant. 

Introducing a dimensionless time 2eVt/h and making 
use of the constraint on the Cfc , one finds that for the cor- 
relation function (H) each term of the perturbation series 
depends on the Josephson coupling only through the 
combination {Ej/V'^~p)'^'"'. Therefore, the perturbative 
expansion in Ej will finally lead to a power series in the 
applied voltage V. This is a consequence of the special 
form of the correlation function at zero temperature. 

It is instructive to first discuss the problem in several 
limits where results are already available. We start with 
the regime of classical phase diffusion which corresponds 
to taking the limit p = 0. Writing the expansion (||) in 
terms of a continued fraction Hi one obtains 
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where Q{x) is the Heaviside step function. This corre- 
sponds to the zero temperature limit of the result ob- 
tained by Ivanchcnko and Zil'berman The current- 
voltage characteristic (0) starts with an ohmic line and 
displays a cusp at = RIc- A further increase of the 
voltage results in a decreasing current. 

As discussed above, an expansion in the Josephson cou- 
pling energy at p = amounts to an expansion in 1/V. 
The series (||) corresponds to the Taylor expansion of (0) 
given by 
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This series has a finite radius of convergence and is re- 
stricted to \V\ > RIc- The limit of convergence just 
coincides with the position of the cusp. Hence, we see 
that for small p and low temperatures the expansion (|^) 
will only converge for sufficiently large voltages. 

To examine this further, we now allow for finite p and 
turn to the regime of charging effects where Ec ^ Ej. 



Within the standard theory of environmental effects on 
Coulomb Blockade (CB) [^, tunneling is treated per- 
turbatively and the current- voltage characteristics is ob- 
tained from the leading term (n — 1) of the series (^) 
For an environment with an ohmic low frequency com- 
ponent, one finds a zero bias anomaly of the Cooper pair 
current Q 
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While for p > 1 this result describes the suppression of 
the current by the Coulomb blockade effect, for p < 1 it 
corresponds to a divergent zero bias conductance. One 
might hope that higher order terms regularize the diver- 
gence, but the discussion of the classical phase diffusion 
limit suggests that for p < 1 the series is divergent for 
small V . We will show below that this is indeed the case. 

Now, for small p and small voltages, corrections to 
the linear part of the current-voltage characteristic (|^) 
arise from Macroscopic Quantum Tunneling (MQT). In 
this regime the Josephson junction is mostly in its zero 
voltage state and the voltage drop occurs at the resis- 
tor. Occasionally, a phase slip will cause a finite voltage 
across the junction leading on average to a finite dc con- 
tribution. In the overdamped limit 2tt^p^Ej <C Ec the 
current-voltage characteristic is given by 
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Note that the zero bias differential conductance goes to 
zero for p < 1 where we found a divergence within CB 
theory. On the other hand, the zero bias differential con- 
ductance of (|l^) diverges for p > 1. 

In order to reconcile these findings we make use of the 
analogy between a Josephson junction and a damped par- 
ticle in a periodic potential. Schmid ||ll| has noted that 
the regions of small and large p are related by a self- 
duality of the model. Further progress [ ^Jl3| in the cal- 
culation of the mobility of the damped particle has been 
based on the thermodynamic Bethe ansatz . Duality 
has also been exploited in the context of the fractional 
quantum Hall effect 12 15|]. 



It should be emphasized that self-duality relies on 
strictly ohmic damping. In contrast, the spectrum (|^) 
has a cutoff frequency lor. This does not spoil dual- 
ity in the long-time limit of the correlation function (|^) 
where the dependence on lor can be absorbed in an ef- 
fective voltage scale introduced below. The long-time 
limit restricts us to small voltages with eV <^ JiujR or 
equivalently V/ RIc ^ Ec/ir^ p^Ej. For the typical case 
of small environmental impedances, voltages of interest 
are of order RIc^ and then a rather wide range of ratios 
Ec/Ej is allowed. Below we will show how the I-V- 
curve can be extended beyond the voltage limit imposed 
by the strictly ohmic approximation. 
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As a consequence of duality, the current voltage char- 
acteristics can be obtained from an integral representa- 
tion 1 13 Iq] which implies two complementary expansions 



describing the zero temperature behavior. In the scaling 
limit the CB series (0) takes the form 
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On the other hand, the MQT series, with the leading 
order terms (Mfl, reads 
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The coefficients of these expansions are determined by 
duality and given by 
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The two expansions (y_l|) and (|12| ) have a finite radius 
of convergence which can be expressed in terms of a crit- 
ical voltage 
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For jO < 1 the series (^T|) converges for V > Vc and thus 
is a large voltage expansion, as discussed above, while 
the series (|l2|) converges for low voltages V < Vc- The 
role of the expansions is interchanged above p — 1 where 
( |rT| ) yields a low voltage expansion while (12) determines 
the large voltage behavior. The situation is illustrated in 
fig. 1^ where on the right side it is indicated which expan- 
sion converges on which side of the curves. 

As we have already seen for the case p = 0, the current- 
voltage characteristic displays a peak. This is still true 
for finite p as seen from fig. ^ where we present /-In- 
curves for various values of p. Here, we show the current 
as a function of the voltage Vj across the Josephson junc- 
tion which is related to the externally applied voltage by 
Vj = V — RI . As can be seen from the figure, for fi- 
nite external impedance, the /-V-curve corresponds to 
a peak of finite width and a maximum current suppressed 
with respect to the critical current in the absence of fluc- 
tuations. As p is decreased, the peak narrows and the 
usual Cooper pair current at zero voltage builds up. 

Due to the complicated form of the expansions ( pT| ) 
and ( p^ ) an analytic determination of the position and 
height of the peak in the current-voltage characteristics 
is in general not possible. However, it turns out that to 
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FIG. 2. Range of convergence of the two expansions ( |l l| ) 
and ^ for [e< /■k'^){Ec/ Ej) = 0.5, 1 and 2 (dotted, full and 
dashed line, respectively). Below the horizontal dashed line 
p = 1, the series (|l^) converges to the right of the curves while 
above p < 1 it converges to their left. The opposite holds for 
the series dl2). 
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FIG. 3. Zero temperature Cooper pair current-voltage 
characteristics for {e< /ii'^){Ec/ Ej) = 1 and p = 0.001, 0.01, 
0.05 and 0.1 shown as full, dotted, dashed, and dashed-dotted 
line, respectively. 

a very good approximation the peak position is given by 
the critical voltage (p^). The quality of this approxima- 
tion can be seen from fig. ^ where we compare it (dashed 
lines) with the peak positions determined numerically 
(full lines). The lines have been restricted to values of 
p for which eVmax < ^-QTvjJr in order to ensure the ap- 
plicability of the theory. Likewise, the maximum current 
can quite reliably be estimated by /max ~ Vc/R- For 
more precise results, a numerical evaluation of (11) and 
( |T^ ) is required which does not present special problems. 

As already emphasized, the range of validity of the ex- 
pansions ([l^) and ( p^ ) is restricted by the assumption 
of strictly ohmic damping. However, the current-voltage 
characteristics can be extended to larger voltages if CB 
theory yields a good description at the limit of validity. 
One may then continue the J-y-curve to larger voltages 
by using the result of CB theory with the full frequency- 
dependent impedance (^. 



3 



1.0 



0.8 



^ 0-6 



0.4 



0.2 



0.0 



0.0 



0.05 



0.1 
P 



0.15 



0.2 



FIG. 4. Comparison of peak positions (full lines) and the 
critical voltage (|l|) for [e'' /%^){Ec/ Ej) = 0.01, 0.1, 1, 10, and 
100 increasing from the upper to the lower curves. Results are 
only shown for values of p which satisfy eVmax < IQTiujr. 



can be determined within the approximation of a self- 
dual model with strictly Ohmic impedance. The results 
of the dual model have been connected with the phe- 
nomena of Coulomb blockade and macroscopic quantum 
tunneling. This allows for the calculation of the I-V- 
curve for experimentally relevant frequency-dependent 
impedances in a large range of parameters accessible by 
state-of-the-art technology. 

The authors would like to thank M. H. Devoret, D. 
Esteve, and A. Steinbach for many inspiring discussions. 
One of us (GLI) is grateful to the SPEC for hospitality 
and the Volkswagen-Stiftung for financial support during 
a stay at the CEA Saclay. The other author (HG) was 
supported by the DAAD through PROCOPE. 



To illustrate this point we consider the case p = 1/2 
which for strictly Ohmic damping allows for an exact so- 
lution 1 17 1^ . Summing up the two expansions (pT| ) and 
(|l^), one finds in both cases 
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Obviously, this result no longer describes a peak struc- 
ture as it was the case for p < 1/2. This is however 
an artifact of the assumption of ohmic damping. For 
p = 1/2 the plateau value of ( |l6| ) just corresponds to 
the value (^) given by CB theory. Beyond the valid- 
ity of the expansions (|ll|) and (^2|) one may thus use 
the cutoff-dependent leading order term of the series (H) 
which describes the decrease of the current for larger volt- 
ages and leads again to a peak in the I-V-cmve. For 
Ec ^ Ej, ( ]l^ ) reaches its plateau value for very small 
voltages. Then, CB theory only fails for voltages below 
{Ej/Ec)RIc- Quite generally, for p < 1 and Ec :S> Ej 
the result of CB theory with the exact impedance can 
be employed except for small voltages. There, however, 
the strictly ohmic approximation is appropriate and the 
results discussed above can be used. 

For p > I, CB theory describes the low voltage behav- 
ior. For example for p = 2, the case dual to p = 1/2, one 
finds 
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which agrees with (|^) to leading order in V. This result 
for strictly Ohmic damping diverges for large V but is in 
fact regularized by the cutoff in the impedance. Thus for 
p > 1 and Ec ^ Ej the entire peak in the I-V-cmve is 
determined by CB theory. 

In conclusion, we have determined the shape of the 
Josephson current peak in the presence of quantum fluc- 
tuation. We have shown that the low voltage behavior 
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